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Abstract 

Winter hardiness is an important trait for grapevine breeders and producers, so identification of the regulatory mechanisms 
involved in cold acclimation is of great potential value. The work presented here involves the identification of two grapevine 
ICE gene homologs, ValCEl and ValCE2, from an extremely cold-tolerant accession of Chinese wild-growing Vitis amurnensis, 
which are phylogenetically related to other plant ICE! genes. These two structurally different ICE proteins contain previously 
reported ICE-specific amino acid motifs, the bHLH-ZIP domain and the S-rich motif. Expression analysis revealed that ValCEl 
is constitutively expressed but affected by cold stress, unlike ValCE2 that shows not such changed expression as a 
consequence of cold treatment. Both genes serve as transcription factors, potentiating the transactivation activities in 
yeasts and the corresponding proteins localized to the nucleus following transient expression in onion epidermal cells. 
Overexpression of either ValCEl or ValCE2 in Arabidopsis increase freezing tolerance in nonacclimated plants. Moreover, we 
show that they result in multiple biochemical changes that were associated with cold acclimation: Vo/C£7/2-overexpressing 
plants had evaluated levels of proline, reduced contents of malondialdehyde (MDA) and decreased levels of electrolyte 
leakage. The expression of downstream cold responsive genes of CBF1, C0R15A, and C0R47 were significantly induced in 
Arabidopsis transgenically overexpressing ValCEl or ValCEl upon cold stress. ValCE2, but not ValCEl overexpression induced 
KIN! expression under cold-acclimation conditions. Our results suggest that ValCEl and ValCE2 act as key regulators at an 
early step in the transcriptional cascade controlling freezing tolerance, and modulate the expression levels of various low- 
temperature associated genes involved in the C-repeat binding factor (CBF) pathway. 
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Introduction 

Grapevine ( Vitis L.) is one of the most widely cultivated fruit 
crops worldwide, and is of great economic importance. However, 
grape production is often severely limited by various biotic and 
abiotic stresses [1-5]. For example, low temperature stress greatly 
restricts the geographic range of grapevine cultivation, and 
decreases berry yield and quality. Significantly, the grapevine 
cultivars that currently dominate the market in terms of acreage 
and production of premium wines are derived from the species 
Vitis vinifera, but they tend to be very sensitive to low winter 
temperatures [6]. Thus, enhancing low temperature tolerance of 
grapevine is of great practical importance. 

In this context, the current study focuses on Vitis amurensis Rupr., 
a wild grape species that is native to China and is extremely cold- 
tolerant [7], withstanding freezing temperature as low as — 40°C 
[8]. V. amurensis therefore has great potential as an experimental 
system to identify mechanisms of cold tolerance and as a 



germplasm resource for grapevine cold-resistance breeding. Some 
highly cold-resistant Vitis cultivars have previously been produced 
using classical breeding methods; however, resistance to cold stress 
is a multigenic trait, which limits the effectiveness of using 
traditional breeding [9]. Therefore, understanding the mecha- 
nisms underlying tolerance and adaptation to cold stress could 
potentially lead to the development of new strategies for improving 
the yield of cold sensitive agronomic plants and expanding the 
geographic areas of production. 

Cold acclimation has been extensively studied, resulting in 
considerable evidences at the molecular level that cold stress 
triggers a multitude of physiological responses [10]. Cold responses 
are complex and highly regulated via activation of signaling 
pathways and numerous genes encoding proteins that act direcdy 
in stress tolerance. To date, numerous cold-regulated (COR) genes 
have been functionally identified, such as antifreeze proteins, later 
embryo abundant (LEA) proteins, molecular chaperones and 
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enzymes involved in detoxification and biosynthesis of osmopro- 
tectants [9], [1 1—13]. Ectopic expression of some COR genes have 
been reported to result in improved cold tolerance, including 
studies with tobacco [14-16], rice [17], [18], strawberry [19] and 
Arabidopsis [20], [21]. Conversely, the expression of Craterostigma 
plantagineum or spinach LEA proteins in tobacco did not induce any 
significant changes in freezing tolerance [15], [22]. Such findings 
are not surprising since a clear increase in freezing tolerance is 
rarely obtained by expressing a single cold-induced gene, even if 
the end product is directly related to development of freezing 
tolerance. 

Recent studies have shown that Arabidopsis cold stress tolerance 
can be enhanced by modulating the signaling pathways triggered 
by low temperature stress [23]— [25] . Such pathways include the 
CBF (C-repeat binding factor), also known as DREB (dehydration 
responsive element binding)//CE (inducer of CBF expression) 
signaling pathway, which have also been characterized in other 
plant species [25]— [2 7] . It was reported that increased expression 
of the entire battery of COR genes resulted from overexpressing the 
Arabidopsis transcriptional activator CBF1 [23], [24], [28], [29]. 
The underlying mechanism involves the CBF1 protein binding to 
the CRT/DRE regulatory element located in the promoters of the 
target genes, thereby inducing cold response metabolic pathways 
and enhancing cold tolerance. Jaglo-Ottosen et al. [23] further 
found that constitutive overexpression of CBF1 induces expression 
of the COR genes in non-acclimated Arabidopsis plants and 
increased freezing tolerance at a whole plant level, an effect that 
was not observed by expressing COR15A alone. This further 
suggests that freezing tolerance is a multigenic trait involving genes 
with additive effects. Overexpression of Arabidopsis CBF genes, 
which reside at the nodes of regulatory networks in cold responses, 
were able to improve chilling/ freezing tolerance in different plant 
species, or homologs from other plant species could enhance the 
freezing tolerance of transgenic Arabidopsis [30], [31]. Other known 
regulatory components include an upstream transcription factor, 
ICE1 (inducer of CBF Expression 1), which encodes a MFC-type 
transcription factor that positively regulates CBF3, and which plays 
a critical role in cold acclimation [32]. Additionally, the Arabidopsis 
gene AtICE2, a positive regulator belonging to the bHLH family, 
has been shown to activate AtCBFl [33]. 

Since the discovery of the Arabidopsis ICE genes, homologs have 
subsequently been found in a variety of crop species, including 
wheat, rice, banana, tea, trifoliate orange and grapevines [34]- 
[44], In addition, 7C£-like proteins have been overexpressed in 
transgenic plants and shown to increase stress tolerance [34], [35], 
[37]. For example, overexpression of an /Ci?-like gene from V. 
amurensis in tobacco was reported to result in a increased cold 
tolerance [42], and Arabidopsis has been shown to become more 
tolerant of cold, drought and salt stresses when expressing the V. 
vinifera ICE la or -lb genes [43]. In this current study, we 
investigated whether two members of the ICE gene family from 
a highly cold-tolerant accession of V. amurensis differ functionally 
from those previously identified from other Vitis species, or 
whether they play similar roles in activating multiple components 
of cold acclimation responses. 

Materials and Methods 

Plant materials, growth conditions, and cold treatment 

Two-year-old Chinese wild-growing Vitis amurensis accession 
'Heilongjiang Seedling' potted plants developing from stem 
cuttings were grown in a greenhouse under natural photoperiod. 
For the cold treatment, plants with a uniform growth status were 
transferred to a chamber (LT-BIX120L, LEAD-Tech (SHANG- 



HAI) SCIENTIFIC INSTRUMENT CO., LTD, China) at 0°C 
with a 16 h photoperiod (200 |J.mol m -2 s _1 light) and 8 h dark. 
Unstressed plants were used as controls (0 h). Leaves from 
untreated (control) and cold-treated V. amurensis plants were 
harvested at time points 0, 1, 3, 6, 12, 24, and 48 h and 
immediately frozen in liquid nitrogen, prior to RNA extraction. 
More than 3 plants were collected and pooled for each time point, 
and the sampling was in triple for biological replicates. 

Structural features, phylogenetic tree and expression 
analysis 

Total RNA derived from grapevine leaves were isolated using 
the plant RNA kit (Omega Bio-tek, Doraville, GA, USA), and first- 
strand cDNA synthesis was synthesized from 2 ug of total RNA 
using PrimeScript RT Reagent Kit according to the manufactur- 
er's manual (TaKaRa, Dalian, China). Two partial lengths of V. 
amurensis ICE cDNA fragments were amplified by PGR using 
degenerate primers 5'-TGGACTSSTCSTCKTCGTGYTCK- 
CC-3' and 5'-TCATCMGCCTCRTCWWTCAAACCSGA-3', 
5 '-CCTYCAGTKGGBKC ACAGCCMACKCT-3 ' and 5'-GC- 
AAACCYATTGAARCAGCTGATRACMGC-3 ' , which were 
designed based on the known nucleotide sequence of ICE 
homologs from other plant species. Isolation of the full length 
cDNA sequences was carried out using the SMART RACE Kit 
(Clontech, Palo Alto, CA, USA). Amplification was performed at 
94°C for 3 min; 27 cycles of 94°C for 30 s, 55 to 58°C for 30 s, 
and 72°C for 1 min; followed by 5 min at 72°C. The cDNA pools 
for 3' and 5' RACE were generated from total RNA extracted 
from cold-stressed leaves of V. amurensis, using the RNAprep Pure 
Plant Kit (Omega). Subsequendy, a nested PCR was performed 
with the prepared cDNA pool using the adaptor primer UPM and 
ValCEl gene-specific primer (GSP) 5'-CATCCA- 
CATGTTCTGTGGCCGCAGACCAGG-3' for 3'RACE, 5'- 
CTGATGACCGCTTGCTGAATGTCTAGCCCAAGGC-3' 
for 5'RACE; and the VaICE2 gene-specific primer: 5'- 
GCCGGCGTCAGACAGTACTGGAAGCTTAGG-3 ' for 
3'RACE, 5'-ATCGCTCGCACTGGAGCTTTGTC- 
GAAGCGC-3' for 5'RACE. Amplicons were cloned into the 
pMD19-T vector (TakaRa) for sequencing and the full-length 
cDNAs of ICE homologs were predicted by comparing and 
aligning the 5'- and 3'-RACE amplified sequences, using BioEdit 
software (Version 7.0.1). The putative full-length cDNAs were 
amplified using primers designed from the extreme 5' and 3' ends 
and 5'RACE-Ready cDNA as template, cloned to pMD-19-T 
(TaKaRa) to generate pMD-ValCEl, and pMD-VaICE2, and 
verified by sequencing. Chromosomal location prediction was 
performed using the BLAT server (http://www.genoscope.cns.fr/ 
cgi-bin/vitis/webBlat) at the Genoscope Genome Browser. The 
molecular weights (MW) and isoelectric points (pi) of the 
corresponding proteins were predicted with the ProParam tool 
(http://www.expasy.ch/tools/pi_tool.html). Nuclear localization 
signals were predicted based on the predicted protein sequence 
using the online server (http://www.predictprotein.org/) and 
homolog searches were conducted with the NCBI BLAST server 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). Molecular model build- 
ing of bHLH-ZIP domain of ValCEl or VaICE2 was carried out 
using SWISS-MODEL server (http://swissmodel.expasy.org). 
Sequence alignment was performed using the DNAMAN software 
(Version 6.0, Lynnon Biosoft) and the phylogenetic tree was 
constructed with MEGA 5.1 software [45] using the neighbor- 
joining method. 

To determine the expression profiles of ValCEl and VaICE2 in 
V. amurensis and Arabidopsis plants, total RNAs were extracted as 
described above. After treatment with RNase-free DNase, the first- 
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strand cDNA was synthesized using the PrimeScript First Strand 
cDNA Synthesis Kit (TaKaRa). Semi-quantitative RT-PCR was 
performed to assess the expression of ValCEl and VaICE2 over a 
cold stress time course. The ValCEl cDNA was amplified using 
the primers 5'-ATGTTACCCAGGTCGAACGACGT-3' and 
5'- CTACAGCATACCGTGGAAGCC TG-3', and the VaICE2 
cDNA was amplified using the primers 5 ' -ATGGTGTC C A- 
GAGTGAACGGCGTC-3' and 5 ' -CTAC AGC ATACC GTG- 
GAAGCCTG-3'. Grapevine GAPDH (GenBank accession 
no. CB973647) was used as a loading control using the following 
primers: 5 '-TTCTCGTTGAGGGCTATTCC A-3 ' and 5'- 
CCACAGACTTCATCGGTGACA-3 ' . Arabidop.su Actin2 
(AT3G 18780) was served as reference control using the following 
primers: 5 ' -CTTGC AC C AAGC AGC ATGAA-3 ' and 5'- 
C C GATC C AG AC ACTGT ACTTC C TT- 3 ' . Three replicates 
were performed for each semi-quantitative RT-PCR reaction. 

Subcellular localization and transactivation activity assay 

To construct green fluorescent protein (GFP) translational 
fusion vectors, full-length ValCEl and VaICE2 cDNAs were PCR- 
amplified with the primers 5 ' -ggaattcCATATG ATGTT AC C - 
CAGGTCGAACGACGT-3 ' (EcoR I) and 5'-gcGTCGACCTA- 
CAGCATACCGTGGAAGCCTG-3 '(&/ I), 5'-catg CCATG- 
GATGCTGTCCAGAGTGAACGGCGTC-3' (Mo I) and 5'- 
cgGGATCCCTACAGCATACCGTGGAAGCCTG-3'(5am// 1), 
and fused in-frame upstream of the of GFP reporter gene in the 
pBI22 1 vector. The empty vector containing only GFP sequence 
was used as a positive control. The isolated plasmids were 
concentrated to ~1 u.g/u.1, and used to coat one set of gold 
particles for bombardment experiments. ValCElv.GFP and 
VaICE2v. GFP fusion proteins were transiently expressed in onion 
(Allium cepa) epidermal cells using particle bombardment and 
subsequent localization of the proteins was performed as 
previously described [46]. Images were collected from the 
transiently transformed onion epidermal cells (n&20) that 
experienced a 24-h dark culture (22°C), using a confocal 
laser-scanning microscope (LAM510, Carl Zeiss GmbH, Jena, 
Germany). 

For the transcriptional activation assay, the sequences corre- 
sponding to the predicted open reading frames (ORFs) of ValCEl, 
amplified with 5 ' -gc TC T AG A ATGTTAC C C AGGTC GAAC - 
GACGT-3' (Xba I) and 5 ' -gc TCTAGA CAGCATACCGTG- 
GAAGCCTGCCG-3' (Xba I), and VaICE2, amplified with 5'- 
ccg C TC G AG ATGC TGTC C AG AGTG AAC GGC GTC - 3 ' (Xho 
I) and 5'-ggGGTACCCAGCATACCGTGGAAGCCTGCCG- 
3'(Kpn I), were fused in frame with the GAL4 DNA binding 
domain in the pGBKT7 vector, to make pGBKT7- ValCEl and 
pGBKT7- VaICE2, respectively. These constructs, as well as pCLl 
(positive control) and pGBKT7 (negative control), were trans- 
formed into yeast strain AH 109 cells and transformants were 
streaked on plates containing SD/Trp- and SD/Trp-/His-/Ade- 
medium. After incubation at 28°C for 3 days, the growth of the 
transformants were evaluated. A /J-galactosidase assay was carried 
out according to the manufacturer's instructions (Clontech). 

Generation Arabidopsis plants overexpressing ValCEl or 
ValCE2 

The cDNA sequences corresponding to the ORFs of ValCEl or 
VaICE2 was amplified from 5'RACE-Ready cDNA using the 
following gene-specific primers containing sequences for restriction 
endonucleases: 5 ' -ccAT_CGAT ATGTT AC C C AGGTC GAAC - 
GACGT-3' (Cla I) and 5'-gc TCTAGA CTACAGCATACCGTG- 
GAAGCCTG-3' (Xba I) for ValCEl and 5'-cc ATCGA- 
TATGCTGTCCAGAGTGAACGGCGTC-3' (Cla I) and 5'- 



ccgCTCGAGCTACAGCATACCGTGGAAGCCTG-3 ' (Xho I) 
for VaICE2. The PCR products were cloned immediately adjacent 
to the 3' end of the CaMV35S promoter in the pART-CAM-S 
vector [47]. The resulting construct was confirmed by sequencing 
and transferred into Agrobacterium tumefaciens strain GV3101 by 
electroporation. Arabidopsis transformation was performed by the 
floral dip method [48]. Tl seeds were collected from individual 
lines and screened on 50 mg/mL kanamycin-contained medium 
to analyze the segregation of the resistant phenotype. T2 
kanamycin-resistant seeds were harvested from non-segregating 
families and confirm the kanamycin-resistant phenotypes. T3 
homozygous lines were validated by RT-PCR using the primers as 
above described in grapevine expression analysis, and further used 
for all experiments. 

Stress tolerance assay 

For the cold treatment, surface-sterilized T3 homozygous seeds 
of vector-carrying control and ValCEl I ^-overexpressing Arabidop- 
sis lines were germinated and grown on MS medium supplement- 
ed with 50 mg/L kanamycin for one week, then transferred to pots 
containing a mixture of perlite: sand: peat (1:1:1, v/v) for two 
weeks. The three-week-old seedlings of the transgenic plants 
overexpressing ValCEl and VaICE2 (n = 50) and control plants 
(n = 50) were transferred to the chamber (LT-BIX120L) at 0°C 
under the condition as above-described in grapevine for 0, 6, 12, 
24, or 48 h. The freezing tolerance assay was performed by 
transferring 3-week-old seedlings (n = 50) grown in pots to the pre- 
chilled chamber (LT-BIX120L) at — 6°C for 8 h under continuous 
dim light conditions (2.5 umol m 2 s ') and subsequendy 
returning them to normal conditions. Survival rates of plants 
were evaluated after 7 days. The freezing treatment experiment 
was performed in triplicate. 

Determination of electrolyte leakage, and 
malondialdehyde and proline content 

Three-week-old seedlings from the empty vector control 
Arabidopsis plants (n = 20), ValCEl and VaICE2 T3 transgenic 
lines (n = 20) were subjected to the cold stress treatment described 
above and leaves were harvested to determine three biochemical 
features associated with cold tolerance: electrolyte leakage was 
assessed by ion leakage analysis as previously described [49]; 
proline content was determined using acid-ninhydrin reagent and 
acetic acid [50]; and malondialdehyde (MDA) content was 
measured using the thio-barbituric acid (TBA) method as 
previously described [51]. 

Analysis of cold-responsive genes regulated by ValCEl or 
ValCEl 

Total RNA was extracted from vector-carrying control and T3 
transgenic Arabidopsis leaves treated at 0°C for 0, 3, 12, 24 and 
48 h using the plant RNA kit (Omega) following the manufactur- 
er's instructions. cDNA was synthesized from 2 u.g of total RNA 
using the PrimeScript RT Reagent Kit with the oligo (dT)18 
primer, according to the manufacturer's instructions (TakaRa) and 
quantitative real-time PCR was performed using the SYBR 
Premix Ex ™ Taqll kit (TakaRa) on an iCycler iQ5 thermal 
cycler (Bio-Rad). The reactions were carried out in triplicate in 96- 
well plates (25 ul/well) in a mixture containing 12.5 ul 2xSYBR 
Premix Ex ™ Taqll, 1 ul each of primer (10 U.M in stock), 1 u.1 
template cDNA and 9.5 ul ddH20. Two-step real-time PCR 
reactions were performed under the following conditions: 95°C for 
10 s, followed by 40 cycles of denaturation at 94°C for 15 s, 
annealing and extension at 57°C for 30 s, and data acquisition at 
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57°C for 15 s. Three replicate PCR amplifications were 
performed for each sample. Transcript levels of AtlCEl, AtICE2, 
AtCBFl, AtCBF2, AtCBF3, AtCOR15A, AtCOR47, AtRD29A 
and AtKINl were measured in ValCEl or VaICE2 overexpress- 
ing or empty vector control plants. The amount of transcript for 
each gene, normalized to the internal reference Atactin2, was 
analyzed using the 2-AACt method [52]. Oligos used for real- 
time PCR were: 5'- C TTC CATC C GTTG AC AC C T AC -3' and 
5'-CTCTAGCTTGCTGGCCTTTAG-3' for AtlCEl; 5'- 
TCCACAAACGCTGTCTTACC-3' and 5'-GTTCACTGCC- 
TTTCCTTCTCT-3 ' for AtICE2; 5 ' -GAGAC GATGGTGGA- 
AGCTATTT-3' and 5'-AGCATGCCTTCAGCCATATTA-3' 
for AtCBFl; 5 ' -GACCTTGGTGGAGGCTATTT-3 ' and 5'- 
ATCCCTTCGGCCATGTTATC-3' for AtCBF2; 5'-GACGT- 
TGGTGGAGGCTATTT-3 ' and 5'-AGCATCCCTTCTGC- 
CATATTAG-3' for AtCBF3; 5'-GGCGTATGTGGAGGAGA- 
AAG-3' and 5'-CCCTACTTTGTGGCATCCTTAG-3' for 
AtC0R15A; 5 ' -GGCTGAGGAGTAC AAGAAC AA-3 ' and 5'- 
ACAATCCACGATCCGTAACC-3' for AtCOR47; 5'-GCAA- 
TGTTCTGCTGGAC AAG-3 ' and 5'- TCCTTCACGAAGT- 
TAACACCTC-3' for AtKINl; 5'-GCTTTCTGGAACAGAG- 
GATGTA-3' and 5'-CGACTCTTCCTCCAACGTTATC-3' 
for AtRD29A. Data analysis shown was obtained from three 
biological replicates. 

Results 

Characterization of two ICE homologs from V. amurensis 

Two /CE-orthologs, designated ValCEl (GenBank accession 
no. KC8 15984) and VaICE2 (GenBank accession no. KC8 15985), 
were isolated from leaves of the cold-tolerant Chinese wild- 
growing V. amurensis accession 'Heilongjiang Seedling'. The 
complete cDNA of ValCEl is 1,949 bp, comprising a 97 bp 
5'UTR, 301 bp 3'UTR and a 1,551 bp ORF that was predicted 
to encode a 516 amino acid protein with a MW of 55.5 kDa and a 
pi of 5.30. The complete cDNA of VaICE2 is 2,108 bp, comprising 
a 214 bp 5'UTR, 277 bp 3'UTR and a 1,617 bp ORF that was 
predicted to encode a 538 amino acid protein with a MW of 
58.38 kDa and a pi of 5.05. A sequence alignment showed that the 
two proteins differ significantly from each other (60% amino acid 
similarity). Chromosome location of the two ICE homologs in the 
V. vinifera cv. Pinot Noir clone P40024 genome suggested that 
ValCEl is mapped on chromosome 1, spanning 2714 bp, while 
VaICE2 is mapped on chromosome 14 and spans 4234 bp 
(Fig. 1A). Analysis of structural properties revealed that the 
predicted ValCEl or VaICE2 protein possesses the typical 
features of ICE proteins, including a serine-rich region (S-rich), 
a basic helix-loop-helix (bHLH) domain, an ICE-specific domain 
[38], a zipper region (ZIP) and an ACT_UUT-ACR-like domain 
(Fig. IB). An interPro scan suggested that both proteins belong to 
the AfTC-like bHLH family of transcription factors, and in support 
of this, they contain nuclear location signals (NLSs) at position 12— 
46 aa for ValCEl and 339-362 aa for VaICE2 (Fig. IB). The 
predicted three dimensional structures of the bHLH-ZIP domain 
of ValCEl and VaICE2 were distinct, and were more similar to 
the structures of DIMER (PDB ID: lr05) and HETERO DIMER 
(PDB ID: 2ql2), respectively (Fig. 1C). 

Phylogenetic analysis of various ICE amino acid sequences 
indicated the existence of two major groups (Fig. 2A). VrICE4 is 
phylogenetically distinct from the other groups in which they could 
be clearly classified into dicot and monocot specific subgroups. 
Our phylogetic data demonstrated ValCEl and VaICE2 fell into 
different clades of the dicot subgroup, with ValCE 1 being closely 
related to the previously reported ValCEl 4 (98.7% identity) from 



V. amurensis [42], VrlCEl (97.7% identity) from V. ripara [44], 
whereas VaICE2 being closely similar to VrICE2 (99.4% identity) 
from V. ripara [44], VvICEl (98.9% identity) and VvICEla (98.7% 
identity) from V. vinifera [43]. In addition, VvICElb shared 
98.72% identity with VrICE3, which belongs to another clade of 
dicot subgroup. Comparison of ValCEl, 2 and their homologs 
AtlCEl and AtICE2 from Arabidopsis, as determined by 49—60% 
sequence similarity, revealed they share highly conserved regions 
in the bHLH DNA binding domain and ACT domain in their C- 
terminal regions (Fig. 2B). In contrast, only a moderate sequence 
conservation was found in the ZIP domains between Vitis and 
Arabidopsis ICE proteins. A potential sumoylation site previously 
identified in Arabidopsis ICE proteins [53] was also observed in 
VaICE2, and the S-rich region, which has been suggested to be a 
site of phosphorylation [32], [54] was present in either ValCEl or 
VaICE2. Finally, the R236 residue that substituted with H236 in 
the Arabidopsis icel mutant [32], causing a loss of ICE function, was 
shown to be present in the two V. amurensis ICE proteins (Fig. 2B). 

To gain insight into the biological functions of ValCEl and 
VaICE2, semi-quantitative RT-PCR was used to examine their 
expression profiles in leaves of V. amurensis over a cold stress time- 
course (Fig. 2C). VaICE2 transcript levels were constant over the 
time course, while a rapid induction of ValCEl expression was 
observed after 1 h, followed by a gradual decline from 3 h to the 
minimum level after 12 h, and finally a return to background 
levels at the end of the time course. These results suggest that 
ValCEl and VaICE2 may be involved in cold stress responses. 

ValCEl /ValCE2 is nuclear-localized and can act as 
transcriptional activators in transient assays 

Sequence analysis showed that ValCEl and VaICE2 have 
putative NLSs, suggesting that they target the mature proteins to 
the nucleus. To test this, ValCEl and VaICE2 were transiendy 
expressed as translational fusions at the N terminus of GFP in 
onion (Allium cepa) epidermal cells. Confocal imaging showed that 
GFP alone (control) was present in the cytoplasm and nucleus, as 
expected, whereas cells transformed with either ValCEl ::GFP or 
VaICE2::GFP showed strong fluorescence exclusively in the nucleus 
(Fig. 3A). 

Yeast one hybrid assays were carried out to determine whether 
ValCEl and VaICE2 possess transactivation activity. As shown in 
Fig. 3B, yeast transformants with pCLl (positive control), 
pGBKT7- ValCEl and pGBKT7- VaICE2 grew well on the SD/ 
Trp- medium, but also grew normally on the SD/Trp-/His-/Ade- 
medium, and exhibited fairly strong /?-galactosidase activity. In 
contrast, transformants carrying the negative control (pGBKT7) 
did not grow on the SD/Trp-/His-/Ade-medium and could not 
show /S-galactosidase activity. Thus, ValCEl and VaICE2 have 
transactivation activities in yeasts. 

Over-expression of ValCEl or ValCE2 in Arabidopsis 
increases cold tolerance 

To determine whether ValCEl or VaICE2 enhances cold stress 
tolerance, Arabidopsis WT-type Col-0 was transformed with 
constructs carrying empty vector, ValCEl- or VaICE2-codmg 
regions under the control of constitutive 35S promoter (Fig. 4A). 
Expression of ValCEl or VaICE2 transgens in these T2 lines was 
confirmed by semi-quantitative RT-PCR, and two ValCEl (L2 
and L3) and two VaICE2 transgenic T3 homozygous lines (LI and 
L6) with similar transcript levels (Fig. 4B) were selected for further 
analyses. These two independent transgenic lines over-expressing 
ValCEl and VaICE2, together with the corresponding empty 
vector-transformed control lines, were subjected to a whole plant 
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Figure 1. Location and structure of the ValCEl, 2 genes. (A) Gene structure and locus. Chromosomal localization of ValCE1,2 genes were 
predicted in V. vinifera cv. Pinot Noir clone P40024 genome. ValCEl comprises four exons spanning 2.7 kb with a predicted 1,551 bp ORF, mapping to 
Chromosome 1 and position 21694886-21700313. ValCE2 contains four exons spanning 4.2 kb and has a predicted 1,617 bp ORF, which is mapped 
to Chromosome 14 and position 24967920-24976387. Nucleotide numbers are indicated above the gene structure. (B) Schematic protein structures 
of ValCEl and ValCE2 with the N-terminal domain containing a S-rich motif and the C-terminal domain represented by helix-loop-helix (HLH), ICE- 
specific domain, Zipper region (ZIP), and ACT_UUR_ACR-like (ACT) domain. The putative nuclear localization signals (NLS) are indicated as green 
boxes. Codon numbers are indicated above the protein structures. (C) The predicted tertiary structures of the bHLH-ZIP domains of the (a) ValCEl and 
(b) ValCE2 proteins evaluated by the SWISS-MODEL sever. 
doi:10.1371/journal.pone.0102303.g001 



freezing assay to evaluate plant survival after freezing treatment 
(Fig. 4C). The freezing test consisted of exposing 3- week-old 
transgenic and control plants to a temperature of — 6°C for 8 h, 
before recovery at normal temperatures. The ValCEl and VaICE2 
transformed plants exhibited less freezing damage and increased 
survival rates compared with empty vector control plants after a 7 
d recovery period (Fig. 4D). Only ~5% of the control plants (LI 
and L2) survived the treatment, while the Fa/C£7-overexpressing 
plants exhibited survival rates of 91% for L2 and 81% for L3, and 
85% for LI and 95% for L6 in Fa/C£ , 2-overxpressing plants 
(Fig. 4D). No obvious differences in phenotype were observed 
between the empty vector control plants and ValCEl or ValCEZ 
transgenic plants grown under normal conditions. These results 
suggest that ValCEl or VaICE2 overexpression leads to enhanced 
freezing tolerance. 

Over-expression of ValCEl or ValCE2 affects electrolyte 
leakage as well as proline and MDA metabolism 

To investigate the physiological and biochemical factors that 
might contribute to the improved cold tolerance of the transgenic 



plants, MDA and proline content were evaluated over a time 
course of exposure to cold stress (0°C), as well as electrolyte 
leakage, which is commonly used as an index of membrane injury 
[23], [55]. A representative experiment comparing the empty 
vector control LI and 3 5S:: ValCEl L2 and 35S::VaICE2 L6 
transgenic plants is shown in Fig. 5. Under non-acclimation 
conditions, electrolyte leakage of all tested lines did not vary 
greatly i.e. between 22% in the control (LI) and 21% in the 
VaICEl-(L2) and VaICE2-(L6) overexpressing lines (Fig. 5A). 
Electrolyte leakage increased with the prolonged cold stress and 
reached 47% in the control and 40% and 39% in ValCEl and 
VaICE2 lines, respectively, after 48 h. 

Proline accumulation in plants has been associated with a wide 
variety of environmental stresses, and confers stress tolerance by 
facilitating osmotic adjustment, protecting proteins and mem- 
branes and quenching reactive oxygen species [56], [57]. 
Significantly, for this study, increases in proline levels have been 
well-documented in Arabidopsis and other plants during cold 
acclimation [58], [59]. Similar changes were observed in both 
ValCEl and VaICE2 overexpressing Arabidopsis lines (Fig. 5B), since 
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Figure 2. Phylogenetic tree, alignment and expression profiles of ValCEl, 2. (A) Phylogenetic tree based on the deduced amino acid 
sequences of ICE from a range of plant species. A neighbor-joining tree, with the following predicted ICE protein sequences, with GenBank accession 
numbers listed in parentheses: V. amurensis (ValCEl, AGP04217; ValCE2, AGP04218; ValCE14, ADY17816), V. vinifera (VvlCEl, AFI49627; VvlCEIa, 
AGQ03810; VvlCEl b, AGQ0381 1), V. riparia (VrlCEl, AGG34704; VrlCE2, AIA58705; VrlCE3, AIA58706; VrlCE4, AIA58707), Arabidopsis thaliana (AtlCEl, 
NP_189309; AtlCE2, NPJ72746), Brassica napus (BnlCEl, AEL33687), Brassica rapa subsp. Chinensis (BrICEl, ACB70963), Capsella bursa-pastoris 
(CblCEl, AAS79350), Camellia sinensis (CslCE, ACT90640), Eucalyptus camaldulensis (EclCEl, ADY68776), Eucalyptus globules (EglCEl, AEF33833), 
Eutrema salsugineum (EslCE, ACT68317), Glycine max (GmlCEl, ACJ39211), Hordeum vulgare (HvlCE2, ABA25896), Malus x domestica (MdbHLHI, 
ABS50251), Oryza sativa (OslCEl, Osl 1g0523700; OslCE2, Os01g0928000), Populus suaveolens (PsICEl, ABF48720), Populus trichocarpa (PtrlCEl, 
ABN58427), Raphanus sativus (RslCEl, ADY68771), Triticum aestivum (TalCE41, ACB69501; TalCE87, ACB69502), Zea mays (ZmlCE2, ACG46593), was 
produced by ClustalX 2.0 alignment followed by tree construction using MEGA 5.0 with 100 bootstrap tests. The branch support values are indicated. 
The length of the scale bar corresponds to 5 substitutions per site. (B) Comparison of ICE amino acid sequences from V. amurensis (ValCEl, 2) and 
Arabidopsis thaliana (AtlCE1,2). Deduced amino acid sequences were aligned using ClustalW. Sequences and accession numbers are shown for the 
following: Vitis amurensis (ValCEl, KC815984; ValCE2, KC815985) and A. thaliana (AtlCEl, AAP14668; AHCE2, AA063441). Residues in black and gray 
regions indicate identical and similar residues, respectively, between isoforms. Four predicted domains are labeled: a S-rich motif, a basic-helix-loop- 
helix-leucine zipper (bHLH-ZIP) region, /Cf-specific domain [38] and a ACT-UUR-ACR-like domain. The red triangles indicate the position of the 
mutation isolated by Chinnusamy et al. [32] and Kanaoka et al. [54], and the residue targeted for sumoylation by SIZ1 [53]. Green triangle indicates E- 
box/N-box specificity site [69], blue asterisks indicate core residues for DNA binding sites [70], and red diamonds indicate dimerization interface/ 
polypeptide binding sites [70], (C) Expression profiles of ValCEl and ValCE2 during a cold stress time-course experiment. Semi-quantitative RT-PCR 
was used to determine ValCEl and ValCE2 transcript levels in cold treated grapevine leaves at indicted times. Grapevine GAPDH was used as a loading 
control. 

doi:1 0.1 371 /journal.pone.01 02303.g002 



similar levels of proline were detected in the transgenic and wild 
type plants without cold stress, but substantial proline accumula- 
tion occurred in cold-treated ValCEl and VaICE2 transgenic lines. 
After 48 h of cold stress, the proline contents in the leaves of 
ValCEl or VaICE2 transgenic lines were ~2. 3-fold higher than 
those of the control plants (Fig. 5B). 

MDA is also considered an indicator of plant oxidative stress 
and structural integrity of the membranes in response to low 



temperature [60], [61]. Under non-acclimation conditions, all 
tested lines had similar MDA levels ranging from 1.76-2. 1 1 ug/ g 
FW (Fig. 5C). However, MDA accumulation was considerably less 
in the ValCEl or VaICE2 transgenic lines than in the controls 
during the cold stress time-course. These results indicated that 
ValCEl and VaICE2 overexpression resulted in increased levels of 
proline, but decreased levels of MDA and reduced levels of 
electrolyte leakage. 
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Figure 3. Nuclear localization and transcriptional activation assay of ValCE1,2. (A) Confocal imaging of the ValCEl and ValCE2 transiently 
expressed in onion epidermal cells as GFP fusion proteins, driven by the 35S promoter. GFP alone was used as a control. Left hand panels show dark 
field images of green fluorescence, middle panels show the morphology of the cells in bright field and the right hand panels show the merged 
images. (B) Fusion proteins of pGBKT7-l/a/CF7, pGBKT7-Va/Cf2, pCLI (positive control) and pGBKT7 (negative control) were expressed in the yeast 
strain AH109. Transformants were incubated on SDATrp- and SDATrp-/His-/Ade- to assess their growth and tested for jS-galactosidase activity. 
doi:10.1371/journal.pone.0102303.g003 



ValCEl or ValCE2 positively regulate cold-induced gene 
expression 

To further identify molecular components associated with the 
ValCEl and VaICE2 mediated stress tolerance that we observed in 
the 35S::VaICEl L2 and 35S::VaICE2 L6 transgenic Arabidopsis 
lines, the expression patterns of genes in the ICE-CBF pathway 
were evaluated by real-time PCR (Fig. 6). Under control 
conditions, AtlCEl or AtICE2 expression was relatively low in 
either the empty vector control or ValCEl- and VaICE2- 
overexpressing plants, respectively. However, overexpression of 
ValCEl or VaICE2 resulted in a reduction in transcript levels 
compared with the control during the cold stress time course. The 
possible effect of ValCEl or VaICE2 on the expression of 
downstream genes through transcriptional regulation was also 
investigated (Fig. 6). Under control conditions, transcripts of none 
of the three CBF genes was detected in the control, but higher or 
similar level CBF transcript levels were present in the ValCEl or 
VaICE2 overexpression lines. Expression of the CBF genes showed 
a general decline in the control plants starting at 12 h of cold 
stress. Starting at 3 h, CBF1 and CBF2 transcript levels were 
considerably greater in the transgenic lines than in the control 
plants, while a substantially greater expression in the transgenic 
lines was not observed for CBF3 until 24 h. Expression of the CBF 
downstream target genes {COR15A, COR47, KLN1 and RD29A) 
were also investigated. COR47 had extremely high expression in 
both transgenic lines and the vector control compared with the 
other genes, and COR15A exhibited higher expression in the 
transgenic lines compared with the control after 48 h. The KTN1 



gene was expressed at considerably higher levels in the VaICE2 
overexpression than that in the ValCEl overexpression line and 
the vector control, while no substantial difference in RD29A 
expression was seen among the three different genotypes. These 
results suggest that ValCEl and ValCE2 positively regulate the 
expression of the CBF genes by differentially controlling down- 
stream genes in response to cold stress, which in turn likely 
contributes to freezing tolerance. 

Discussion 

In the current study, two putative /CE-orfhologs, ValCEl and 
VaICE2, from the highly cold-tolerant species V. amurensis, were 
structurally analyzed and functionally tested. We report that these 
two genes encode MYC-type bHLH transcription factors that are 
nuclear-localized, and are candidate regulators of the CBF gene 
expression during cold stress. Reports on /CB-orthologs in three 
genotypes of Vitis species [42]-[44] might indicate that there are 
multiple /C£-like genes in grapevine. However, sequence com- 
parisons revealed that none of the mRNA sequences that were 
previously reported [42], [43] completely matched the two 
sequences identified here (Fig. SI in File SI). Thus, Vitis ICEs 
were categorized on the basis of sequence similarity (Fig. S2A in 
File SI), phylogenetic clustering (Fig. S2B in File SI) and genetic 
locus (Table SI in File SI). Our results provide evidence of the 
presence of 4 polymorphic loci for the identified ICE-like genes in 
grapevine. ValCEl, ValCEl 4 and VrlCEl are alleles from which 
they are located in grapevine chromosome 1, whereas VaICE2, 
VrICE2, VvICEl and VvICEla are alleles, mapping to chromosome 
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Figure 4. Freezing tolerance evaluation of 35S::ValCE1,2 transgenic Arabidopsis plants. (A) A schematic map of the T-DNA region of 
35S::ValCE1 ,2 fusion constructs employed for Arabidopsis transformation. RB, right border; LB, left borders; CaMV35S, Cauliflower mosaic virus 35S 
promoter; OCS, octopine synthetase terminator; NOS, nopaline synthase promoter; NPTII, Kanamycin resistance gene; NOS-T, nopaline synthase 
terminator. (B) RT-PCR was used to assess the transcript abundance of ValCEl or ValCE2 in the transgenic Arabidopsis plants. Atactin2 was used as a 
reference control (C) 3-week-old plants were treated at — 6 C C for 8 h and then transferred back to normal conditions for recovery. Photographs were 
taken after 7 d of recovery. (D) Survival rates of plants exposed to — 6°C. Average survival rates and standard errors were calculated using the results 
of three separate experiments with 50 seedlings per line for each freezing stress. 
doi:10.1371/journal.pone.0102303.g004 




Figure 5. Effect of ValCEl, 2 expression in Arabidopsis on levels of (A) electrolyte leakage, (B) proline and (C) malondialdehyde 
(MDA). Three-week-old Arabidopsis plants of vector-carrying control, ValCEl L2, and ValCE2 L6 were grown at 0°C for the time indicated. Leaves were 
collected to assess electrolyte leakage and free proline and MDA. Each value is the mean ± SD of three replicates in case of electrolyte leakage, 
proline and MDA content (6 seedlings each) for the indicated time points. 
doi:10.1371/journal.pone.0102303.g005 
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Figure 6. Effect of ValCE 7,2overexpression in Arabiclopsis on transcript levels of genes involved in the cold stress pathway. qRT-PCR 
analysis using leaves from control Arabidopsis LI, 35S::ValCE1 L2 and 35S::ValCE2 L6 plants. Three-week-old plants were grown at 0 r C for the time 
indicated. The tested genes were /If/CFJ, /4f/C£2, AtCBFI, AtCBF2, AtCBF3, AtCOR15A, AtCOR47, ATKIN1 and AtRD29A. Atactin2 was used as a reference 
control. Each value is the mean of three replicates. 
doi:10.1371/journal.pone.0102303.g006 



14. Additionally, the other independent cluster that VvICElb and 
VrICE3 involved are alleles, which are located in chromosome 17. 
Of these three distinct categories, very few, and mostly conserved, 
substitutions were observed for each category, but it appears that 
VrICE4 is phylogenetically distinct from the above-mentioned 
three categories and located in chromosome 18. A detailed 
sequence assay on VrICE4 revealed that it lacked the typical S-rich 
motif of ICE proteins, presented 3 amino acid mutations (A-T, A- 
T, and E-D) in ICE-specific domain, and exhibited remarkable 
sequence difference in ACT_UUR_ACR-like (ACT) domain (Fig. 
S2A in File SI). Moreover, result from FLAGdb""" database 
prediction indicated that VrICE4 has the best match with 



GSVIVG0 1009234001 with functional annotation as DNA 
binding protein rather than inducer of CBF expression 2 (Table 
SI in File SI) when compared with other grapevine ICE proteins. 
Therefore, further investigation on VrICE4 is still required to 
confirm its role as /CE-like transcription factor. 

Evidence from structural and phylogenetic analyses further 
suggested that ValCEl and VaICE2 may have different properties 
when involved in cold stress. This is also suggested by their 
expression profiles (Fig. 2C), which showing that they are both 
constitutively expressed, but that the expression of ValCEl is 
affected by cold treatments. With the availability of the reported 
expression data on grapevine ICE genes [42-44] , a comparison of 
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cold responses of these genes on different genotypes revealed that 
ValCl, 2, VvICEla, b and VrICEl-4 were all detected before and 
after cold treatment, which is most consistendy observed to be 
constitutively expressed in other plants [32], [34], [40]. However, 
VaICE14 is a notable exception where the transcript was not 
detected under non-stress condition but under low temperature 
[43]. One possible explanation for such difference has been 
proposed [44] , in which this discrepancy might be attributed to the 
different genotypes and sampling materials derived from different 
cultural ways. An additional effort in validation of the expression 
profile of this gene before and after cold stress in repeated trials is 
required to clarify its expression property. 

Many of the physiological and biochemical changes that occur 
during cold acclimation are direcdy correlated to an up-regulation 
of COR gene expression, which is activated by the constitutive 
expression of the ICE transcription factor. Liu et al. [62] found 
that overexpression of AtlCEl in cucumber was sufficient to 
increase chilling tolerance and to simultaneously alter the levels of 
several cold responses associated factors (e.g., free proline, MDA 
and soluble sugars). Li et al. [43] discovered that the survival rates 
of VvICEla and VvICEl A-overexpressing Arabidop.sis lines were 
significandy higher than those of the wild type under cold stress. 
Similarly, overexpression of Solarium lycopersicum ICE1 was reported 
to improve tolerance of chilling stress in tomato plants, as 
indicated by differences in electrolyte leakage [37]. In our study, 
overexpression of either ValCEl or VaICE2 in transgenic 
Arabidopsis plants increased freezing tolerance at the whole-plant 
level, based on survival rates (Fig. 4B, C), which correlated with 
increased accumulation of proline, and a reduction in MDA and 
electrolyte leakage (Fig. 5A-C). These results suggest that although 
an endogenous ICE-CBF pathway is present in Arabidopsis system 
and that enables the ectopic VaICEl,2 transgene expression to 
influence the plants' freezing-tolerance capacity. There are also 
considerable evidences to suggest that proline [63]-[65], MDA 
[14], [66], [67], and electrolyte leakage [25], [68], contribute to an 
enhancement of freezing tolerance. Our results confirmed that the 
greater tolerance of ValCEl or Fa/GEP-overexpressing plants to 
cold stress positively linked to elevated proline levels, and 
negatively correlated to the reduction of MDA content as well as 
electrolyte leakage. However, these data did not reveal a difference 
between the ValCEl and VaICE2 transgenic lines. It is important 
to note that plants overexpressing ValCEl or 2 had similar levels of 
proline and MDA contents relative to the vector-carry controls 
under normal growth conditions (0 h), but upregulated in proline 
level and downregulated in MDA content under time-course cold 
stresses. We therefore assume that the commonly observed 
accumulation of free proline and MDA under non-acclimation 
are partially determined by proline biosynthesis gene (e.g. P5CS, d- 
l-pyrroline-5-carboxylate synthetase) or MDA reductase gene 
itself property with a relatively low expression, but also induced by 
the expression of a transcription factor under cold-stressed 
conditions. Additionally, ionic leakage from the cells is considered 
as an indicator that the semipermeable nature of the plasma 
membrane has been lost, at least transiendy, in response to 
freezing. This is probably the main cause for the similar electrolyte 
leakage observed in controls and ValCEl-/ Fa/CE?-overexpressing 
lines under non-stressed conditions. 

We propose that our analyses of the ValCEl or VaICE2 
overexpressing Arabidopsis plants do not contradict the general cold 
acclimation model [32], [43], but rather suggest a more integrated 
circuitry involved in cold acclimation signaling. Chinnusamy et al. 
[32] reported that over-expression of AtlCEl increased the 
expression of AtCBF2, AtCBF3 and cold-regulated genes (CORs) 
under cold stress. In contrast, Fursova et al. [33] found that over- 



expression oiAt!CE2 only resulted in dominant changes in AtCBFl 
transcription levels after cold acclimation, while Badawi et al. [34] 
showed that over-expression of wheat ICE genes in Arabidopsis 
induced a higher expression of AtCBF2, AtCBF3 and some COR 
genes only after cold acclimation. Our data suggest that the 
expression of AtlCEl or AtICE2 was slightly lower in the ValCEl 
and VaICE2 overexpressing Arabidopsis plants. One possible 
explanation for this is that co-suppression of the endogenous gene 
occurred in the target plant, but we suggest that the V. amurensis 
and Arabidiopsis ICE1 or ICE2 cDNA nucleotide sequences are not 
highly homologous (64.6-69.9%). The observation that several 
downstream genes were significandy induced in ValCEl- or 
Fa/C£2-overexpressing lines might suggest that the constitutive 
expression of the transgenes affects the expression of the 
endogenous gene. Thus, as ValCEl and VaICE2 mRNA levels 
are abundant in the transgenic plants, competition at a post- 
transcriptional level could explain this slight decrease in endog- 
enous AtlCEl or AtICE2 transcript levels. We speculate that 
another signal transduction pathway may exist in the CBFl 'DREB1 
gene regulatory network that is stimulated by ValCEl or VaICE2. 
It should be noted that the activation of the expression of these 
downstream genes differed between the ValCEl and VaICE2 
overexpressing lines (Fig. 6). The enhanced tolerance of 
35S::VaICEl and 35S::VaICE2 plants coincides with an up- 
regulation of the stress-responsive genes CBFl, COR15A and 
COR47. However, some differences were seen in the expression 
patterns between the transgenic lines, and particular in the case of 
KLN1, which suggests functional differences between ValCEl and 
VaICE2. In the case of VvICEla and VvICEl b transgenic Arabidopsis 
lines, both contribute to the modulation of AtRD29A and AtCOR47 
in response to cold stress [43] . Together with the above-mentioned 
examples, these differences in targeting to the downstream genes 
are likely due to a varying number of ICE homologs in different 
plant species [32], [33], [40-44]. However, it has not yet been 
established which of the identified grapevine ICE genes specifically 
control the different sets of cold-responsive genes. A detailed 
comparison of the expression of these ICE genes and their target 
genes through microarray analysis might help address this issue. 

Taken together, ValCEl and VaICE2 are two previously 
unreported ICE-like transcription factors from V. amurensis, whose 
regulatory roles in cold acclimation are suggested by the results 
presented here. Both genes may act as positive regulators to 
increase the levels of cold-responsive genes in transgenic lines 
under cold stress. Moreover, ValCEl and VaICE2 influence cold 
stress-related factors such as electrolyte leakage, and levels of 
proline and MDA, thereby alleviating damage by ROS and 
enhancing osmotic protection. Our data may help elucidate the 
cold-acclimation pathways of Vitis species and more ultimately 
guide the design of strategies for improving the stress tolerance of 
agricultural crops. 

Supporting Information 

File SI Contains the files: Figure SI Multiple alignment 
of the mRNA sequences of /C£-homologous from 
different Vitis species. Identical nucleotide sequences are 
highlighted on a black background while white boxes indicate at 
least three identical nucleotides. The GenBank accession numbers 
are reported as follows: ValCEl (KC8 15984), VaICE2 
(KC815985), ValCEH (HM231151), VvICEl (JQJ707298), VvI- 
CEla (KC831748), VvICElb (KC831749), VrlCEl (KF994961), 
VrICE2 (KF994962), VrICE3 (KF994963) and VrICE4 
(KF994964). Figure S2 Protein sequence similarity and 
phylogenetic clustering of ICE from three different Vitis 
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species. (A) Protein sequence alignment of 10 grapevine ICEs. 
Identical residues are outlined in black. Amino acids are 
numbered on the right. (B) Phylogenetic tree based on the 
deduced amino acid sequences of ICEs from three different Vitis 
genotypes. A maximum-likelihood phylogenetic tree of the amino 
acid sequences of ICE from V. amurensis (ValCEl, AGP04217; 
VaICE2, AGP04218; VaICE14, ADY17816), V. vinifera (VvICEl, 
AFI49627; VvICEla, AGQ03810; VvICElb, AGQ0381 1), and V. 
nparia (VrlCEl, AGG34704; VrICE2, AIA58705; VrICE3, 
AIA58706; VrICE4, AIA58707) is constructed by MEGA 5.0 
with 1000 bootstrap tests. The branch support values are 
indicated. The length of the scale bar corresponds to 0.05 
substitutions per site. Table S 1 Feature lists of BLAST 
results or queries of ten grapevine ICE genes available 
in FLAGdb ++ . A Blastp search with an E- Value of l.E-50 was 
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